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1.  INTRODUCTION 


At  the  Harry  Diamond  Laboratories  (HDL)  magnetic  signature  facility  at  Aberdeen,  MD,  magnetic 
signatures  of  tanks  are  determined  by  measuring  the  field  in  a plane  under  the  tank  at  points  spaced  1 ft 
apart  in  a 20  x 40  ft  square  grid  (J.l.  Cooperman  and  T.  H.  Hopp,  unpublished  report).  Forty-one  readings 
are  taken  (one  for  each  foot)  as  the  tank  advances  over  a linear  array  of  21  magnetometers.  The  problem 
treated  here  is  how  best  to  correlate  different  sets  of  measurements  to  eliminate  statistical  iineertainties  in 
the  measurement.  In  a given  run,  the  tank  may  take  a path  that  is  slightly  shifted  with  respect  to  its  pre- 
vious run,  and  the  process  of  recording  readings  on  the  magnetometers  for  every  foot  of  travel  of  the  tank 
across  the  linear  array  may  have  begun  at  a slightly  different  time.  One  cannot,  therefore,  simply  compare 
the  output  of  corresponding  magnetometers  and  treat  the  problem  by  the  usual  statistical  methods.  There 
are  good  reasons,  owing  to  the  small  shifts  described  above,  that  corresponding  magnetometers 
give  the  same  reading  during  different  runs. 

A second  problem  related  to  these  measurements  is  sorting  out  the  tank’s  distortion  of  the  earth’s 
magnetic  field  from  the  larger  residual  field  of  the  tank.  The  latter  is  easily  measured  by  bucking  out  the 
earth’s  field  with  large  Helmholtz  eoils  existing  at  the  site,  but  the  former  must  be  inferred  by  comparing 
these  measurements  with  similar  measurements  when  the  earth’s  field  is  not  bucked  out.  Shifts  similar  to 
those  described  above  must  again  be  considered  when  one  compares  these  measurements. 

Both  problems  are  treated  here  in  a unified  manner  by  considering  the  question,  “What  can  be  said 
about  a function  of  two  variables  if  it  is  sampled  only  at  a finite  number  of  points  in  a square  grid  covering 
a given  area?”  This  problem  is  solved  by  expanding  the  matrix  of  measurements  in  a discrete  Fourier  series 
with  the  largest  number  of  low-frequency  components  that  can  uniquely  be  determined  from  the  measure- 
ments. This  expansion  provides  a natural  and  uncontrived  extension  of  the  function  continuously  and 
smoothly  to  any  point  within  the  area.  By  comparing  the  functions  generated  in  this  manner  from  two 
separate  runs,  one  may  then  continuously  shift  the  functions  with  respect  to  each  other  to  obtain  a best 
fit.  In  this  way,  an  analytic  determination  of  the  previously  unknown  shifts  between  the  two  runs  is  ob- 
tained, and  systematic  errors  are  eliminated  that  result  when  such  shifts  are  ignored.  Assuming  the  shifts  to 
be  small  compared  with  the  1-ft  spacings  between  data  points,  a closed  expression  for  the  shifts  is  obtained. 
Larger  shifts  may  be  determined  by  using  well-known  computer  methods  for  minimizing  the  explicit  ex- 
pression derived  for  the  square  deviation  between  the  functions. 

The  second  problem  may  be  treated  similarly  by  comparing  the  two  runs  with  and  without  the  earth’s 
field  bucked  out,  except  there  is  a question  as  to  how  the  comparison  should  be  made.  One  can  only 
assume  that  the  correlation  between  the  tank’s  residual  field  and  its  distortion  of  the  earth’s  field  is  small 
(if  not  zero)  compared  with  its  self  correlation,  so  also  a best  fit  between  the  shifted  functions  representing 
the  two  runs  should  be  required  in  this  case.*  The  remaining  field  after  the  best  fit  is  subtracted  is  then  the 
distorted  earth's  field,  provided  it  may  be  assumed  also  that  the  magnetic  behavior  of  the  tank  is  linear  for 
field  strengths  considered  in  these  experiments.  Although  the  formalism  developed  in  this  report  may  be 
useful  in  considering  possible  nonlinear  effects,  th^s  problem  is  not  treated  here. 

In  formulating  the  mathematics  to  handle  the  problem  discussed  above,  the  larger  question  has  been 
answered  of  determining  a function  of  two  variables  by  sampling  at  a discrete  number  of  points  in  a regular 
grid.  The  methods  in  this  report  therefore  more  generally  apply. 

*in  a private  communication,  T.  //.  IIopp  of  UDL  suggests  that  there  is  reason  to  believe  the  distorted  earth's  field 
has  a component  that  correlates  in  phase  with  the  tank's  residual  field.  This  correlation  provides  even  stronger  fustificathn 
for  requiring  a best-fit  nuitch. 
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2.  DISCRETE  FOURIER  EXPANSION 


Let  a sequence  of  readings  on  N detectors  taken  at  M equal  time  or  space  intervals  be  stored  in  an 
NXM  array  A,  where  A(pq)  is  the  q'^  reading  of  the  p^^  detector.  For  convenience,  start  the  numbering  at 
zero  so  that 


0<p<N- 1 , 

(la) 

0<q<M - 1 . 

(lb) 

In  this  problem,  A(uv)  is  a function  of  two  variables,  u and  v,  known  only  at  integer  values  of  its  arguments 
over  the  finite  ranges  given  above.  To  approximate  the  function  for  arbitrary  uv,  A(pq)  is  expanded  in 
finite,  discrete  Fourier  series  and  then  generalized. 

Define  the  Fourier  components  in  the  expansion  of  A by 

N-l  M -1 

A(pq)  = X X e e 

n = 0 m = 0 

(2) 

where,  from 


P=o 


^N5. 


(3) 


expression  (2)  may  be  inverted  to  give 


(NM) 


-1 


N-  1 M-1 

X X 

p=0  q = 0 


.27T 

-1  N pn 


.27T 

-‘m  qm 


A(pq)  . 


(4) 


To  be  consistent  with  this  expansion,  it  is  assumed  that  for  p and  q outside  the  ranges  given  by  expressions 
(la)  and  (lb),  A is  defined  by 


A(p  + N,  q)  = A(pq)  = A(p,  q + M)  . (5) 

Henceforth,  it  is  assumed  also  that  the  implied  summations  in  Dp  and  Sq,  respectively,  extend  over  the 
ranges  given  by  expressions  (la)  and  (lb)  and  similarly  in  and  unless  otherwise  noted.  Now  define 
the  extension  of  A to  arbitrary  uv  by  substituting  p = u and  q = v in  expression  (2)  and  using  the  definition 
of  a^^^  given  by  expression  (4).  Thus, 


—1  V f f NT  - P)*i 

A(uv)  = (NM)  ^ z-  A(p  q ) Z-  e ^ e 

p q nm 


^ X A(p'q')  (p',  - u)  (q;  - v) 


t t 

p q 


(6) 
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where 


Cv,  (p.z)=N“'sin  7rz-^-5r 

sin-^(P  + z) 


= N*"*  sin  7TZ 


cot  “ (p  + z)  + 1 


(7a) 

(7b) 


The  functions  are  a convenient  form  for  evaluating  the  summations  needed  later.  In  apj^ndix  A,  sums 
over  products  of  cotangents  are  discussed,  and  results  are  given  from  which  useful  properties  of  the  can 
be  derived.  These  properties  of  the  Cj^  are  given  in  appendix  B. 


From  expression  (6),  it  is  clear  that  A(uv)  is  not  necessarily  real  for  noninteger  u and  v.  Therefore, 
define  R(uv),  the  extention  of  A(pq)  to  noninteger  arguments,  as  the  real  part  of  expression  (6),  to  give 


R(uv)  = (NM)  * sin  ttu  sin  ttv  ^ A(p^q^) 

/ t 

P q 


7T  / ^ . 77  , / , 

Sin  — (p  --  u)  sin  (q  v) 


(8) 


where  this  expression  has  been  obtained  by  substituting  expression  (7a)  for  the  . Letting  u = p + x/2, 
V = q + y/2,  the  following  expression  results: 


R(p  + x/2,  q + y/2)  = (NM)  ^ sin-^  sin A(p*q') 


[\  p'  p \/2  y q'  q y/2l 

cos  7Tn  + + — + 

, r,  [2  N 2 M J 


p q 


77  , 

Sin  — (p 
N 


• (9) 


x/2)  sin  •—  (q 
M 


>72) 


This  form  is  convenient  for  examining  the  relationship  of  R to  A for  small  deviations  x,y.  As  x and  y ap- 
proach zero,  R(pq)  A(pq). 


Similarly,  if  S is  the  analogous  extension  of  a second  set  of  measurements  stored  in  an  array  B(pq), 
then  the  two  extensions  R and  S may  be  compared  by  shifting  S in  the  opposite  direction  and  minimizing 
their  square  deviation.  In  this  way,  the  two  sets  of  measurements  A and  B are  on  equal  footing.  Thus, 
define 

Q(x,y)  = (NM)'*  ^ rR(p  + x/2,q  + y/2)-S(p-x/2,q-y/2)1  (10) 

pq  L J 

where  S is  given  by  letting  x and  y go  to  -x  and-y  and  substituting  B for  A in  expression  (^).  Carrying  out 
the  algebra  according  to  the  cotangent  sums  in  appendix  A (or  more  simply,  the  Cj^  relationships  given  in 
appendix  B),  one  gets 


CKx,y)  = (1  + cos  77(x  + y)]Qj+  [cos  Try  - cos  7r(x  + y)] 

+ [cos  77X-  cos  7r(x  + y)]  Q3  + [1  + cos  7r(x  + y)”  cos  77x-  cos  Try] 
^ [1  + cos  77(x  + y)j  sin  ttx  sin  Try  Q^(x,y) 

“ [ 1 - cos  TTX  + sin  77X  sin  7r(x  + y)j  sin  Try  Q^(y) 

' [ 1 - cos  Try  + sin  Try  sin  7r(x  + y)j  sin  ttx  Q^(x) 

- [( 1 - cos  TTxX  1 " cos  Try)  - sin  ttx  sin  Try  cos  7t(x  + y)j  Qg 
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where 


Q,  =-^(nm)’‘  I [A2(pq)+B2(pq)j,  (12a) 

Q2='5N'^M'*  rA(pq)A(p'q)  + B(pq)B(p'q)  1 , (12b) 

pp  q L J 

X,  f'A(pq)A(pq')  + B(pq)B(pq')'|  , (12c) 

^ pqq  L J 

Q4=^(a2  + b2),  (12d) 

Q«  (x,y)  = (NM)  ' X cot  ^ (p  - x)  cot  ^ (q  - y)  P(pq)  , (12e) 

Q,(y)=(NM)-‘  1 cotj(q-y)  IP(pq) , (12Q 

Q7(x)=(NM)'*  X cot-^(p-x)  XP(pq).  (12g) 

Qg=AB,  (12h) 

and  where 

A=(NM)’’  I A(pq),  (13) 

pq 

and  the  correlation  matrix  P is  defined  by 

P(pq)=(NM)'*  A(p'q')  B(p' + p,  q' + q)  . (14) 

P q 


In  principle,  the  problem  is  now  straightforward.  Expression  (1 1)  must  be  minimized  with  respect  to  x and 
y— say  this  occurs  at  x = Xo,  y = yo*»  ^hen  the  most  probable  function  F (given  the  two  sets  of  measure- 
ments) is 

F(uv)  = i [^R(u  + xj2,  V + yj2)  + S (u  - xj2,  v - y^/2)]  (15) 

with  an  rms  deviation  a given  by 


a = 


Q(Xo.y 


.)] 


1/2 


(16) 


Minimizing  Q(xy)  may  be  a formidable  task  even  for  a computer  if  N and  M are  of  appreciable  size;  so 
the  problem,  although  well  defined,  is  not  necessarily  satisfactorily  solved  at  this  stage.  In  addition,  there 
may  be  several  sets  of  measurements,  in  which  case  the  procedure  outlined  above  must  be  performed  several 
times  and  the  results  somehow  correlated.  In  the  following  section,  the  case  is  treated  where  Xq  and  y^  may 
be  assumed  small  in  comparison  with  unity.  Then  the  quantities  of  interest  may  all  be  expanded  for  small 
xy,  resulting  in  a considerable  simplification  of  the  problem. 
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Expressions  for  Xq  and  are  given  for  obtaining  these  quantities  by  direct  calculation,  and  no  com- 
puter minimization  procedure  is  required. 

By  use  of  this  procedure,  the  problem  of  treating  more  than  two  sets  of  measurements  is  no  longer 
formidable.  A method  for  correlating  several  sets  of  measurements  to  obtain  the  most  probable  function 
F(uv)  is  described  in  section  4.  The  method  treats  equally  each  set  of  measurements  and  generates  a set  of 
standard  deviations  for  the  shifts  by  which  one  may  estimate  the  errors  in  the  measurements  and  the  degree 
to  which  the  resultant  function  F is  meaningful. 

3.  RESULTS  FOR  SMALL  DEVIATIONS 

If  the  deviations  between  two  sets  of  measurements  are  small,  then  one  may  expand  the  results  of  the 
previous  section  for  small  x,y  and  analytically  minimize  Q(xy).  This  may  be  a good  preliminary  calculation, 
in  any  event,  even  if  the  deviations  are  not  known  to  be  small  or  cannot  be  assumed  to  be  small,  since  the 
resultant  x^,  y^  may  be  used  as  starting  values  to  minimize  the  full  expression  for  Q given  by  expression 


(11). 


Expanding  expression  (1 1)  to  second  order  in  x and  y,  one  obtains 


Q(xy)=  X 

nm 


(17) 


where 


koo  = (NM)  ^ V A(pq)-B(pq) 


pq 


(18a) 


' Z cot  P(pO)  , 

P=1  N 


(18b) 


M-1 


koi=-2M 


I cot^P(Oq), 
q = l ^ 


(18c) 


U8d) 


P(Oq) 


M-1 


(18e) 
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(180 


+ + ^ P(pO)-M-'I  P(0q) 


-2(NM) 


N-l  M-l 

■'  S I 

p=l  q=l 


TTP  7rq  V 

cot  jq-  cot  - P(pq)  , 


and  where  Q2  and  are  given  by  expressions  (12b)  and  (12c)  and  P(pq)  is  defined  by  expression  (14). 

Differentiating  expression  (17)  with  respect  to  x and  y and  separately  setting  the  results  equal  to  zero 
give  the  solutions 

ikoi  - 2kioko2)  / (^*^20*^02  *"  1) 

y0  = 7T^  ^ki  ikjo  - 2kojk2oj  / (^*^20*^02  “ 1)  * (1^^) 

Thus,  by  direct  calculation,  we  have  determined  the  Xq,  y^-that  minimize  Q.  If  the  resultant  values  are  less 

than  0.16  (ttx  < 0.5),  the  sine  and  cosine  expansions  used  to  obtain  these  results  approximate  the  true 

values  to  within  1 percent,  so  the  values  for  Xq,  y^  calculated  by  this  method  may  be  assumed  accurate.  In 
the  limiting  case  in  which  A and  B are  equal,  it  can  be  shown  that  k^^  and  k^ ^ vanish,  so  x and  y are  zero 
as  required. 

The  expansion  for  R for  small  x,  y is 

R(p  + x/2,q  + y/2)^  A(pq) 


-(S)  ^ A(pq')coti(q'-q) 

"Ho*- '> s is] 

^ A(pq')  j(N-l)  ^ + [m  - 1 - cot2  ^(q'  - q)]  ^ 


Try 

2M 


^ TTX  Try 
2N  2M 


A(p'q')  |cot-^  (p'  - p)  cot  ^ (q'  - q)  - l| 


(20) 


where  S means  the  p'  = p is  excluded  from  the  sum;  the  same  holds  for  q.  This  expression  may  not  be 
much  simpler  to  evaluate  than  the  rigorous  expression  (9),  unless  Xq  and  y^  are  sufficiently  small  so  the 
quadratic  terms  may  be  ignored.  In  this  case, 
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F(pq) - T 


A(pq)+  B(pq) 


^ B(p'q)j  cot  ^ (p'  p) 

T r , ,1  rt  , 

— j^A(pq  ) B(pq  )J  cot  - (q  q) 


The  rms  deviation  between  the  shitted  functions  from  expressions  (16)  and  (17)  is  given  by 

1/2 


(21) 


(22) 


CORRELATION  OF  SEVERAL  MEASUREMENT  SETS 
Suppose  tlicre  are  K sets  of  measurements  Aj(pq)  where 

Ki<K 


(23) 


and  each  set  is  shifted  by  some  unknown  amount  from  its  presumably  true  set  of  values  according  to  the 
discussion  in  the  previous  sections.  Let  the  vector  quantity  r-j  represent  the  pair  of  values  where 

these  arc  the  shifts  calculated  to  minimize  the  difference  between  Aj  and  Aj.  That  is,  the  minimum  square 

(24) 


deviation  between  A^  and  A^  obtains  for  the  extended  version  of  Aj,Rj,  evaluated  at 


and  for  the  extended  version  of  Aj,Rj,  evaluated  at 


Then,  clearly 


ated  at 

Rj  (P  - x,j/2,  q - y../2)  . 


(25) 


Now,  suppose  the  centroid  of  all  the  measurements  is  designated  by  the  subscript  zero,  and  for  ex- 
ample, is  the  shift  that  must  be  given  to  Aj  to  bring  it  into  line  with  the  centroid.  Then 


is  the  vector  criterion  for  the  centroid.  To  solve  for  a given  rj^.  it  is  clear  that,  ideally. 


f.o  = '.j  + 'jO  • 


Thus  expression  (27)  becomes 


jo 


V 

i^j 


iO 


V 


= Kr,^  + ^ r.;  = 0 


i^j 


(27) 


(28) 
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or 


(29) 


by  use  of  expression  (26).  These  Tjq  satisfy  expression  (27),  even  though  expression  (28)  may  not  be 
strictly  correct.  It  must  be  stressed  that  expression  (28)  is  rigorously  correct  only  if  the  shifted  matrices 
represent  the  true  function  that  would  be  measured  at  the  shifted  position  and  no  errors  existed  in  the 
measurements.  This  assumption,  of  course,  is  only  approximately  correct. 

Similarly,  although  each  calculated  r^  ideally  should  equal  r^^  + r^^  for  any  n,  taking  r^.  = 0 for  all  i, 
errors  in  the  measurements  and  in  the  extended  functions  prevent  this  equality  from  being  exact.  Devia- 
tions from  equality  therefore  give  some  measure  of  these  errors.  Accordingly,  define  a.j(x)  by 

ai.(x)=|(K-2)'^  I + (30) 

and  define  Oy(y)  similarly. 

Practice  will  determine  how  useful  the  a.j  are,  but  it  is  clear  that  if  ajj(x)  is  not  small  compared  with 
|x-j|,  then  either  some  sel  of  measurements  that  contributed  most  to  the  magnitude  of  ajj(x)  should  be  cast 
out  or  the  whole  procedure  of  shifting  sets  of  data  is  questionable  for  the  given  data.  If  the  a.j  are  reason- 
ably small,  however,  then 

F(uv)=k'  Z Rj(u  + Xj0>v  + yj(j)  (31) 

is  the  best  extension  of  the  data  where  the  rj^  are  given  by  expression  (29).  For  the  case  of  two  sets  of 
measurements  (K  = 2),  expression  (31)  reduces  to  expression  (15)  derived  in  section  2. 
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APPENDIX  A.- COTANGENT  SUMS 


From  E.T.  Whittaker  and  G.N.  Watson,  A Course  of  Modern  Analysis,  (4th  ed.,  Cambridge  University 
Press,  1958,  p.  140,  sect.  7.7,  example  1), 


one  may  sum  over  0 < p < N — 1 to  give 

n cot  ^(p  - a.)  = -N  cot  - a,)  cot  tta^,  + N cos 


ttK  (A*1) 


N 


(A-2) 


where  one  assumes  a ^ a for  n m and  no  a is  an  integer.  If  some  of  the  a are  integers,  the  limits  of 
n m m ^ m ^ 

expression  (A-2)  as  each  approaches  its  integer  value  give 

ricot^(p  -a;)=-N  X cot^(a„,-a.)cottta,„ 


V 

p 


+ X X II  cot -(p  a)fi+cot^— /p  — a.^+Ncos—  (A-3) 

where  extends  over  0 < p < N — 1 excluding  integers  mod^^  /a  ) for  integer  where  mod^ 

(V  R V/ 

a^j  for  n m,  where  v extends  over  the  noninteger  a^^,  and  where  v extends  over  the  integer  a^^. 

‘ ' m ni 

From  expressions  (A-2)  and  (A-3),  then,  for  K = 1 , 

X (P  “ ^)  “ ^ ' ( A-4a) 

Pi) 

and  for  K = 2,  where  a^  ^ a2  integer,  pj  ^ p2» 

cot  ^^p  “ cot  ^ _ 32^  = — N N cot  ~ ^aj  — a-,^  jeot  rra^  - cot  rra^J  , (A-5a) 

p^p  (p  Pl)  ‘'"'n  (p  "2)=  N + Ncot-J  (pi-32)  cot7ra2  + l+cot-^  (p,  a,)  ,(A-5b) 


p^p-^p  ^“'^(P  - Pl)  ^(P  - P2)  = -N  + 2 [1  + ^(P,  - P2)]  ■ (A-5c) 
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Differentiating  expression  (A-4a)  with  respect  to  “a,”  one  may  generate 

K 


y 


cot 


r(p-) 


fK<*) 


where  f satisfies 


with 


^K+l  ^K-l 


fQ  - N;  fj  = — N cot  Tra 


and  f ^ = df/da.  Thus,  solving  expression  (A -7)  for  K = 1 , 


A 

p 


cot^(p-a) 

2 2 
= -N  + N^ 

2 

1 + cot  Tra 

_ _ 

_ 

for  example.  Limiting  expression  (A-9a)  as  a pj  gives 


I 

P^P  ] 


Thus,  the  generalization  of  expression  (A-5c)  is 


^ / 

\ 

cot-(l 

’-P]) 

_ 

_ 

^=3(n-ixn-2). 


I 

p^p 

P=^P 


^ cot  ^ (p  - P,)cotI(p  - p^)=  -N  + 2 1 + cot'^  (p,  - P2) 

2 

+ ^(n2+2)5  . 

3'  ' P 1^2 


(A-6) 


(A-7) 


(A-8) 


(A-9a) 


(A-9b) 


(A- 10) 


Expression  (A-6)  may  be  generalized  by  differentiating  expression  (A-2)  times  with  respect  to  a^ 
for  each  n to  develop  a sum  of  the  form 


in 

p m 


/ 

‘’m 

cot—  (l 

) “ a \ 

= f 

/a. 

. . . 

NV 

mj 

Ql  . 

cr 

By  use  of 

d 7T  / \ 7T 

•^cot— (p  — a)=— 
da  N\  ) N 

the  following  recursion  expression  for  the  generalized  f is  obtained. 


1 . 2 7T  / \ 

1 + cot  -(p  - a) 


(A-11) 


(A- 12) 


= -f„ 


X N d f 

qi-  . .q„  + 1 . . . q^...q^-l  . . . q._  + 


■r  * * 


k Trq^  da„ 
n 


(A- 13) 
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where  the  boundary  conditions  for  f if  each  q = 0 or  1 are  obtained  from  expression  (A-2). 
procedure  is  taken  to  evaluate  expression  (A-11)  if  some  of  the  a^^  are  integers.  The  result 
a^^  = integer,  say  = p^^,  is  obtained  by 


z n 

P’^Pn  IV, 


Similarly,  additional  limits  may  be  taken  if  several  of  the  a^^  are  integers.  Note  that  if  p^^>  N 
a^^>  mod^^Pj^^  , since  ^ ^ (^n)  i^^^P^f^dently  for  each  a^^. 


For  K = 2, 


where 


and 


Z 

P 


/ 

7T  / \ 

cot—  (r 

3 — a,\ 

cot-(p-a2) 

NV 

V 

N\  V 

foo=N;fio  = -Ncot7raj  , 


fl,=-N-Ncot^(aj-a2) 


cot  Tra^  — cot  7ra2 


f20=-N  + N^ 


1 + cot"^  Tra 


= _Nr.ot  Tra^cot^^  a^) 


+ N 


cot  Tra^  ^ "n  (^1  ^2) 


7T 


f22  ~ ^2) 


cot  TTa^  cot  7ra2 


. 2 7T 
sin  VT 


n("i  -"2) 


2 .^,2  7T 


+ N cot 


^ 2 2 
2 + cot  Traj  + cot  7ra2 


A limiting 
for  a given 

(A- 14) 
t . then  let 

(A- 15) 

(A- 16) 

(A- 17a) 
(A- 17b) 

(A-17c) 
(A- 17(1) 

(A-17e) 
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APPENDIX  B,  - THE  (p.z)  FUNCTION 


In  the  form 


Cj^(p,z)  = N * sin  TTZ 


cot  - (p  + z)  + i 


(B-1) 


the  following  properties  of  these  functions  may  be  derived  from  a limiting  procedure  of  he  cotangent 
sums  of  appendix  A: 


Cn(p.0)  = 5p0,  (B-2) 

Cn(p',P  + z)  = (-)'’Cn(p'  + P,z).  (B-3) 

cJj(p.z)  = Cn(-P. -z)  (B-4) 

= Cfj(p,z)- 2N  ■'  isinrrz  , (B-5) 

X Cj^(P.z)=  X C^(-P,  z)=  X Cfj(P  + P^z)  = e”^^  (B-6) 

P P P 

X C^{p-  P',  z)  c*  (p  - p",z')  = (p"  - p'.  z - z')  , (B-7) 


X (p  - p',  z^Cj^  (p  - p",  z'j=  Cfj  (p"  - p',  z - z')  + 2iN’  ' e'^^  sin  rrz' 
= — p” , z'  — + 2iN  * sin  irz  . 
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